Quantum simulator for the Schwinger effect with atoms in bi-chromatic optical lattices 
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Ultra-cold atoms in specifically designed optical lattices can be used to mimic the many-particle 
Hamiltonian describing electrons and positrons in an external electric field. This facilitates the exper- 
imental simulation of (so far unobserved) fundamental quantum phenomena such as the Schwinger 
eflect, i.e., spontaneous electron-positron pair creation out of the vacuum by a strong electric field. 
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Introduction There are several fundamental predic- 
tions of relativistic quantum field theory which have 
so far eluded a direct experimental verification. One 
prominent example is the Schwinger [1] effect (histori- 
cally more accurate would be the name Sauter-Schwinger 
effect, see [2] and [3]), i.e., the spontaneous creation of 
electron/positron- pairs out of the vacuum by a strong 
electric field. For a constant electric field i?, the leading- 
order e^e~ pair creation probability scales as [1~3] 
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where £'s = M c^/{hq) is the critical field strength de- 
termined by the elementary charge q and the mass M of 
an electron (or positron). The above expression (1) for 
Pe+e- does not permit a Taylor expansion in q, i.e., it 
is inherently non-perturbative and thus cannot be repre- 
sented by any finite set of Feynman diagrams. 

Unfortunately, our theoretical understanding of this 
non-pcrturbative QED effect is still very incomplete. 
Apart from the constant field case, only very simple field 
configurations where the electric field either depends on 
time E(t) or on one spatial coordinate such as E(x) are 
fully solved. For example, recently it has been found 
that the occurrence of two different frequency scales in 
a time-dependent field E{t) can induce drastic changes 
in the (momentum dependent) pair creation probability 
[4, 5]. Moreover, the impact of interactions between the 
electron and the positron of the created pair, as well as 
between them and other electrons/positrons is not un- 
derstood. This ignorance is unsatisfactory not only from 
a theory point of view but also in view of planned exper- 
iments which envisage field strengths not too far below 
the critical field strength Eg and thus could be able to 
probe this effect experimentally [6]. 

These considerations motivate the investigation of the 
Schwinger effect via a different line of approach. By suit- 
ably designing a laboratory system, we could reproduce 
the quantum many-particle Hamiltonian describing elec- 
trons and positrons in an electric field and thereby obtain 
a quantum simulator for the Schwinger effect. This would 
facilitate the investigation of space-time dependent elec- 
tric fields such as E{t,x) and should also provide some 
insight into the role of interactions. It should be stressed 



here that our proposal goes beyond the simulation of the 
(classical or first-quantized) Dirac equation on the single- 
particle level, see, e.g., [8-14], but aims at the full quan- 
tum many-particle Hamiltonian. A correct description 
of many-body effects such as particle-hole creation (in- 
cluding the impact of interactions) requires creation and 
annihilation operators in second quantization. There are 
some proposals for the second-quantized Dirac Hamilto- 
nian [15-20] but they consider scenarios which are more 
involved than the set-up discussed here and aim at differ- 
ent models and effects. Similarly, the recent observation 
of Klein tunnelling in graphene [21, 22] deals with mass- 
less Dirac particles - but the mass gap is crucial for the 
non-perturbative Schwinger effect, cf. Eq. (1). 

The model We start with the Dirac equation [7] de- 
scribing electrons and positrons propagating in an elec- 
tromagnetic vector potential which are described by 
the spinor wave- function 5" {h ~ c = 1) 



7^(ia^-g^^)«'-A/«' = 0. 



(2) 



For simplicity, we consider 1+1 dimensions (p, = 0, 1) 
where the Dirac matrices satisfying = 2ri^'^ 

can be represented by the Pauli matrices 7*^ = (T3 and 
7^ = — i(Ti. Since in one spatial dimension there is no 
magnetic field we can choose the gauge qAo = <!> and 
Ai = 0. As a result, the Dirac equation simplifies to 



(3) 



In one spatial dimension, there is also no spin, hence the 
wave function has only two components 5* = (^^,5"^). 
The Hamiltonian for the classical Dirac field reads then 



H 



dx 'iJ^-ia2dx + M<y?. + '^)^ ■ 



(4) 



As the next step, we discrctize the space dimension and 
introduce a regular grid (lattice) Xn = na with a positive 
grid (lattice) constant a and integers n G Z. The dis- 
cretization of the wave function \l/„(i) := -^/a ^'(i, a;„), de- 
fined now at the grid points Xn , gives rise to a discretized 
derivative V^c'^«'(<, x„) [*n+i(t) - *„-i(i)]/(2a). 
Finally, replacing the x-integral by a sum, we obtain 



H 



(5) 
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In order to obtain the quantum many-body Hamiltonian, 
we quantize the discretized Dirac field operators via the 
fermionic anti-commutation relations 

{nit),i'i{t)} = 0, {*,^(i),[*f„W]t} = <5„™<5"^ (6) 

Using = (a„,6„), i.e., 4'^ = a„ and = 5„, the 
discretized many-particle Hamiltonian reads 



H = 



„a„+i + h.c. 



+ 



+ E + ^^)4an + (^n - M)bibn\ • (7) 

n 

The first term describes jumping between the neighbor- 
ing grid points while the remaining two terms can be 
treated as a combination of external potentials. Due to 
the specific form of the jumping, the lattice splits into 
two disconnected sub-lattices: (A) containing a2„ and 
62T1+1 and (B) containing a,2n+i and b2n with integers n. 
Since the two sub-lattices behave basically in the same 
way, it is sufficient to consider only one of them, say A. 
Identifying C2„ = (-l)"a2„ and C2„+i = (-l)"+^62n+i, 
we obtain the form of the well-known Fermi-Hubbard 
Hamiltonian for a one-dimensional lattice 



H =-J 



(8) 



with hopping rate J = l/(2a) and on-site potentials 
V2n = $2n + M and V2n+i = <i>2ri+i - ■ Tlus Hamilto- 
nian will be the starting point for the design of the optical 
lattice analogy. But before we proceed, we note that the 
free part Hq of this Hamiltonian, i.e., with $„ = 0, can 
be explicitly diagonalized. Performing a discrete Fourier 
transform on the lattice 



for p G (— 7r/2a, 7r/2a) and introducing new operators 



1 



we obtain the diagonalized form 



Ha = / dpE, 



ili„ - bIb. 



with the energy spectrum 



Ep = \j AP + ^ cos2(ap) . 



(10) 



(11) 



(12) 



It approximates the relativistic energy-momentum rela- 
tion at the edge of the Brillouin zone, for p ss ±7r/(2a). 
The spectrum of consists thus of two symmetric in- 
tervals separated by a gap of 2M. In order to obtain a 



positive Hamiltonian, we perform the usual re-definition 
Wp O Bp which corresponds to changing the vacuum 
state by filling all Bp states by a fermion. This is com- 
pletely analogous to the Dirac sea picture in quantum 
electrodynamics. In terms of this analogy, A^^ or Ap cre- 
ate or annihilate an electron whereas Bp or create 
or annihilate a hole in the Dirac sea - which is then a 
positron. An additional potential $„, if sufficiently local- 
ized in space, will not modify this spectrum but may in- 
troduce isolated eigenvalues with eigenstates correspond- 
ing to bound states localized in space. 

Experimental set-up The Fermi-Hubbard Hamilto- 
nian (8) can be realized with ultra-cold fermionic atoms 
in a one-dimensional optical lattice with the potential 



W{x) = Wo sin2(2fcx) 4- Al^sin^(fcx) , 



(13) 



where k = -k/ (2a), by taking Wq ^ ISW . Similar settings 
have already been obtained experimentally [25]. 
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FIG. 1: Sketch of the dispersion relation (in units of J and 
1/a) for ISW — (dashed blue curve) and small AW > 
(solid and dotted red curves). 

Although discretization of W{x) at Xn = na gives di- 
rectly Vn of (8), the small perturbation AW doubles the 
original periodicity of the potential W{x) and thus some 
mathematical caution in treating AW perturbativcly is 
required. By a version of WKB approximation for peri- 
odic potentials [24] , we re-derive the Hubbard model and 
the energy band structure from first principles, using two 
sets of Wannier functions localized in the "upper" and the 
"lower" minima of the potential W{x), respectively. An 
interesting universal phenomenon occurs, see Fig. 1. For 
AW = 0, the lowest band is described, in the nearest- 
neighbor approximation used here, by Ep = — Jcos(ap) 
(the dashed blue curve in Fig. 1). Switching on a small 
AW > 0; which immediately doubles the period of the 
potential, forces the energy dispersion relation to halve 
its period (the Brillouin zone shrinks by a factor of two) 
while keeping similar functional dependence on p when 
AW <^ Wq. It leads to splitting of this band into two 
sub-bands (red solid and dotted curves in Fig. 1) which 
are approximately described by 



E±{p) = ±y/AP + J2 cos2(ap) 
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This nearest-neighbor approximation reproduces the re- 
lation (12) for J = l/(2a) and holds uniformly for all 
values of the quasi-momentum p as long as /S.W and J 
are small [26]. In the vicinity of the minimum of the 
upper band and the maximum of the lower band, sepa- 
rated by a gap 2M « ISW, it reproduces the relativistic 
energy-momentum relation. The corresponding Hamil- 
tonian has the same form as the one for the discretized 
Dirac equation (11), thus completing the analogy. 

Using again the WKB approximation, we find that 
the hopping rate J is mainly determined by the poten- 
tial strength Wq, the laser wave- number k and the mass 
TOatom of the atoms (not to be confused with the effective 
mass M of the Dirac particles to be simulated) 

J«.lv^exp|-^yg|, (14) 

where Er = fc2/(2matora) = 7i'^/(8matoma^) is the recoil 
energy. The correction /S.W generates the mass gap M k, 
ISW jl. The analogue of the e^e~ pair creation can be 
simulated if the involved scales obey the hierarchy 

Wosc > J > A/ > T . (15) 

First, the local oscillator frequency cJqsc in the potential 
minima must be larger than J to ensure the applicabil- 
ity of the single-band Fermi-Hubbard Hamiltonian (8). 
Second, J ^ M is required for the continuum limit, i.e., 
that the discretized expression (5) provides a good ap- 
proximation. Similarly, the change A$„ = — of 
the analogue of the electrostatic potential $„ from one 
site to the next should be much smaller than M. Over 
many sites, however, this change can well exceed the mass 
gap 2M, which is basically one of the conditions for the 
Schwinger effect to occur. Finally, the effective temper- 
ature T should be well below the mass gap 2M in order 
to avoid thermal excitations. 

Let us insert some example parameters. For ^Li atoms 
in an optical lattice made of light with a wavelength of 
500 nm, the recoil energy Er would be around 7 fiK. If 
we choose the potential strength as Wq = 10 /iK, the 
hopping rate J would be around 5 /iK which is still suffi- 
ciently below the local oscillator frequency cjosc of around 
34 //K. With AW = 1 //K we would get a mass M of 
500 nK and the effective temperature should be below 
that value. 

Bose-Fermi mapping Since it is typically easier to 
cool down bosonic than fermionic atoms, let us discuss 
an alternative realization based on bosons in an optical 
lattice. To this end, we start with the Bose-Hubbard 
Hamiltonian which has the same form as the Fermi- 
Hubbard Hamiltonian (8) after replacing the fermionic 
Cn by bosonic d„ operators, but with an additional on- 
site repulsion term U {dl^dn ~ l)dl^dn- For large U ^ J 
(which can be controlled by an external magnetic field 



via a Feshbach resonance, for example), we obtain the 
bosonic analogue of "Pauli blocking", i.e., at most one 
particle can occupy each site. Neglecting all states with 
double or higher occupancy, we can map these bosons 
exactly onto fermions in one spatial dimension via 

Cn = dn Y\. exp (^-in dl^drn^ . (16) 

m<n 

As a result, we obtain the same physics as described by 
the Fermi-Hubbard Hamiltonian (8). 
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FIG. 2: Sketch (not to scale) of the temporal stages of the sim- 
ulation (from top to bottom). The solid red curves represent 
the optical lattice potentials W{x) as a function of position x 
and the dashed blue curves correspond to the effective electric 
potential $. The blue solid dots are particles and the empty 
(blue) circle is a hole. 

Simulation procedure The above established analogy 
between the (discretized) quantum many-particle Hamil- 
tonian of electrons and positrons in an electric field, on 
the one hand, and the (Bose or Fermi) Hubbard model 
describing ultra-cold atoms in optical lattices, on the 
other hand, enables laboratory simulations of relativis- 
tic phenomena of strong-field QED. Probably the most 
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prominent of them is the spontaneous pair creation in 
strong electric fields known also as the Schwinger effect 
[1-3] which has to date not been confirmed experimen- 
tally. The original Schwinger effect [1] deals with a con- 
stant electric field E and would correspond to a static 
tilted optical lattice with $(a;) = Ex. However, in view 
of the planned experiments [6], an electric field which is 
localized in space and time is more realistic. In order 
to clearly distinguish non-perturbative spontaneous pair 
creation (via tunneling) from other perturbative effects 
such as dynamical pair creation, the electric field should 
be slowly varying (compared to the frequency scale of the 
gap 2M). 

As a result, we envisage an experimental realization 
sketched in Fig. 2. In order to prepare the initial state, we 
start with AW ^ Wq where the two bands are separated 
by a large gap. The Dirac sea then corresponds to the 
state where all lower minima are filled with atoms while 
all upper minima arc empty (half filling of the lattice, see 
first picture in Fig. 2). If we then decrease AW adiabat- 
ically until AW Wq and thus J ^ il/, the atoms be- 
come delocalized - but still the lower band is filled while 
the upper band remains empty (second picture in Fig. 2). 
As the next stage, we slowly switch on an additional po- 
tential $ to facilitate tunneling from the lower band to 
the upper band - the analogue of the Schwinger effect 
(third picture in Fig. 2). Finally, after slowly switching 
off the potential <I> again (fourth picture in Fig. 2), we 
increase AW adiabatically until AW ^ Wq. By this 
energetic separation, a particle-hole pair created via tun- 
neling is transformed to an atom in one of the upper min- 
ima and, consequently, a missing atom (i.e., hole) in one 
of the lower minima (fifth picture in Fig. 2). This could 
be detected by site-resolvcd imaging [23], for example. 

Note that the creation of a particle-hole pair separated 
by An lattice sites requires the simultaneous tunnelling 
of An particles since two particles are not allowed to 
occupy the same site. This again emphasizes the many- 
particle character of our proposal which goes beyond the 
simulation of the classical Dirac equation. Apart from 
investigating the pair creation probability for space-time 
dependent electric fields E{t,x), this quantum simula- 
tor for the Schwinger effect could provide some insight 
into the impact of interactions. For example, including 
dipolar interactions of the atoms, we may even switch 
between attractive and repulsive interactions by aligning 
the atoms spins parallel or perpendicular to the lattice. 
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